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Abstract 

Scattering of neutrons in the 24-150 keV incident energy range from H2O rel- 
ative to that of D2O and H2O-D2O mixtures was reported very recently. Studying 
time-of-flight integrated intensities, the applied experimental procedure appears to 
be transparent and may open up a novel class of neutron experiments regarding the 
"anomalous" scattering from protons, firstly observed in our experiment at ISIS in 
the 5-100 eV range. The keV-results were analyzed within standard theory, also 
including (1) multiple scattering and (2) the strong incident-energy dependence of 
the neutron-proton cross section aniEo) in this energy range. The analysis reveals 
a striking anomalous ratio of scattering intensity of H2O relative to that of D2O of 
about 20%, thus being in surprisingly good agreement with the earlier results of 
the original experiment at ISIS. 

Keywords: neutron Compton scattering, attosecond physics, quantum entangle- 
ment 



1 Introduction 

Several neutron Compton scattering (NCS) experiments on liquid and solid sam- 
ples containing protons or deuterons show a striking anomaly, i.e. a shortfall in the 
intensity of energetic neutrons scattered by the protons; cf. tD|2|[3]l4|. E.g., neu- 
trons colliding with water for just 100 — 500 attoseconds (1 as = 10~^* s) will see 
a ratio of hydrogen to oxygen of roughly 1.5 to 1, instead of 2 to 1 corresponding 
to the chemical formula H2O. Due to the large energy and momentum transfers 
applied, the duration of a neutron-proton scattering event (the so-called scattering 
time) is a fraction of a femtosecond which is extremely short compared to the usual 
condensed-matter relaxation times. 

Recently this new effect has been independently confirmed by electron-proton 
Compton scattering (ECS) from a solid polymer fT,'4','5^|. The similarity of ECS 
and NCS results is striking because the two projectiles interact with protons via 
fundamentally different forces, i.e. the electromagnetic and strong forces. 

Due to its novelty and far-reaching consequences, however, this effect has been 
the focus of various criticisms, cf. (6l|7l. Therefore, considerable work to identify 
possible sources of experimental and data analysis errors was made during the 
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last five years, which succeeded to demonstrate the excellent working conditions 
of the spectrometer Vesuvio; see |8 | for an account in detail. Extending these 
investigations, the complete "exact formalism" of data analysis (6j was applied 
to NCS-data by Senesi et al. |9|, for the first time. Analysis of time-of-flight 
(TOP) spectra from solid HCl revealed the existence of an "anomalous" decrease 
of the scattering intensity from H by 34%. Moreover, this result was found to 
be in excellent agreement with the corresponding outcome of the standard data 
analysis procedure applied at ISIS (g)- Additionally, the mentioned standard NCS- 
data analysis method | 8 | was successfully compared with a newly proposed (by B. 
Dorner, ILL) model-free data-analysis procedure, the latter being independent of 
the form of the momentum distribution and the resolution function |10| . 

2 NCS in the keV range 

In a novel experiment by Moreh et al. 11 IJ, scattering results of neutrons from 
H2O relative to that of D2O and H2O-D2O mixtures were reported, in the inci- 
dent energy range about 24-150 keV. This experiment was carried out to search 
for the aforementioned anomalous decrease |T| in the neutron scattering intensity 
from protons (relative to that from deuterons) at interaction times in the attosecond 
regime. It is important to note that the energy range here is about 3000 times larger 
than that of our original neutron Compton scattering (NCS) experiment carried out 
with the electron volt spectrometer eVS (newly: Vesuvio) at the neutron spallation 
source ISIS 1 1 1. 

In clear contrast to the original NCS results (T|, Moreh et al. claimed that 
the results in keV range exhibited no anomalous behavior. It was concluded that 
within an overall statistical accuracy of 3% there is no evidence for any deviation 
from the ratios of scattering intensities conventionally calculated on the basis of 
the tabulated total neutron cross sections 1 1 1 1. In addition, it was correctly argued 
that one would have to shake some well established notions in physics to explain 
the aforementioned scattering anomaly. 

Our theoretical analysis of these experimental data, and in particular their com- 
parison with predictions of standard NCS theory 1 12|, is considered in the follow- 
ing. This analysis reveals a striking anomalous ratio of scattering intensity of H2O 
relative to that of D2O of about 20%. Extending the first analysis of single scat- 
tering events 1131 , it will be shown that (a) neither multiple scattering (b) nor 
the strong dependence of the proton total cross section aniEo) on incident en- 
ergy Eo, in the considered keV range, do considerably affect the anomaly under 
investigation. 

2.1 Experimental 

First of all, it is important to notice that both setups, i.e. that of the new keV- 
neutron experiment at the Rensselaer Polytechnic Institute (RPI) 1 1 1 1 and that of 
Vesuvio at ISIS 1 8 1, are basically similar and thus the interpretation of their results 
ought to be based on the same basic theory. The following related remarks should 
be helpful: 

(j) Both are so-called inverse geometry time-of-flight (TOP) setups; i.e., the 
final energy Ef of the measured neutrons is fixed, the neutron initial energy Eq 
is "continuous", and the scattered neutrons are analyzed using a "filter" (of Pe at 
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RPI, with Ef = 24.3(±1.1) keV) or "analyzer foil" (of Au at ISIS, with Ej- = 
4.91(±0.14) eV). 

(ii) The range of scattering angles 9 is similar in both setups, (i.e. one detector 
integrating over 9 = 25° —65° at RPI; 32 detectors measuring at various scattering 
angles in the range 35° - 67° at ISIS). 

(Hi) As a consequence, and according to standard theory 1 12], the characteris- 
tic neutron-proton scattering time in the keV range is shorter by a factor of about 
70 nil as compared to that of the Vesuvio-setup. In view of some of the men- 
tioned theoretical models (see below), however, this difference of scattering times 
may not preclude the appearance of the considered effect in the keV range. 

(iv) The Impulse Approximation (lA) can be safely assumed to be exact in the 
keV range fTll . and it is already known to be sufficiently fulfilled in the eV range 
of Vesuvio l8l ll2l . In simple terms, each neutron scatters from a single nucleus (of 
H, D or O). 

(v) The setup at RPI cannot resolve the neutrons scattered from different nuclei 
and yields instead an integrated scattering signal arising from the (O, H and/or 
D) nuclei of the liquid samples. The ISIS setup, however, provides a H-recoil 
peak well resolved from that of D or O, but this does not represent any significant 
difference that would prevent comparison of results. (E.g., one can simply integrate 
over the individual peak intensities, when required). 

(vi) Furthermore, note that the scattering intensity, measured by a single de- 
tector, is represented by the area of the 24.3 keV line in the TOF spectra, i.e. by 
adding up the total number of counts in the 24.3 keV peak 1 1 1 1. An independent 
fission detector, placed inside a separate beam tube, was employed as a neutron 
flux monitor and served to normalize the TOF spectra. After subtracting the back- 
ground from each signal, the intensity ratios are taken, as presented in Fig. 3 of 
[111. The measurements were repeated about five times, always giving about the 
same results within statistics (personal communication). 

(vii) Already here it should be emphasized that, as already stressed in Ref. tTTl . 
multiple scattering effects do not affect the ratio of scattering intensities from the 
two samples; see also below. 

2.2 Single scattering events 

Recognizing the crucial importance of these novel keV-experimental results, we 
analyzed in detail [13| the data processing scheme indicated in Ref. IlII on the 
basis of standard scattering theory at large energy transfers |8, I2|, where the lA 
is valid; cf. point (iv) above. Some omissions in [ 1 1 J have been revealed and then 
amended, thus leading to a radical revision of the main finding and conclusion of 
Ref. fill. In accordance with the preceding point (vii), in that calculation f 131 
we considered single scattering events only. That calculation 1 13], which also took 
into account the strong dependence of the proton total cross section (th{Eo) on 
incident energy Eq IiT4I . contains no fitting parameter. 

A full account in detail of these investigations is presented in Ref. |13| . The 
results are summarized in Fig. 1 and lead to a radical revision of the main finding 
and conclusion of Ref. |11|: The correct keV-data reduction reveals a strongly 
anomalous ratio of scattering intensity of H2O with that of D2O of about 20%, 
thus being in surprisingly good agreement with the associated results of the original 
ISIS-experiment [JJ. 
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Our results f 131 were confirmed independently by Monte Carlo calculations of 
Mayers [15], applying the routines for NCS-data analysis available at ISIS. Also 
these calculations contain no fitting parameters. 




Figure 1: Intensities of single scattering events ifTJl . Measured (full squares with error 
bars) — taken from Ref. fTT| — and corrected calculation (solid line, open circles) 
for 9 = 45° of scattered intensity ratios versus Xjj, the D2O concentration in the 
H2O-D2O mixture. The discrepancy between experimental data and prediction of con- 
ventional theory is clearly discernible for data with > 0.5. The calculated intensity 
ratio for pure H2O and pure D2O is ca. 23% higher than the experimental one. 



2.3 Double scattering events - Energy dependence of an 

The thickness of all samples in the keV experiment was d = 0.18 cm. However, as 
already mentioned under point (vii) above, is was shown that multiple scattering 
effects did not affect the considered ratio of scattering intensities; explicitly: "The 
effect on the ratio of scattering intensities from the two samples is <I% and was 
neglected."; see page 4 of Ref. |l 1| . 

Applying the Monte Carlo routines available at ISIS, J. Mayers analyzed the 
multiple scattering effects on water and D2O, for the sample thickness of d = 0.18 
cm |I5| . In these calculations, however, the proton total cross section an was kept 
constant (e.g. 15.7 barn). The obtained results clearly showed that, in the H2O 
case, the single scattering events are attenuated due to the considerable sample 
thickness. (Note that the transmittance of the H2O sample (for Eo — 48.6 keV) 
is ca. 81%, see numerical calculations below). But, the crucial point here was 
that the twice scattered neutrons were found to go mainly at forward scattering 
angles, and that 6 « 45° is at the center of the the double scattering events in the 
forward angle range. Thus, as sample attenuation increases, multiple scattering 
also increases. As a result, and including multiply scattered neutrons, the ratio of 
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scattering intensity of H2O relative to tiiat of D2O was calculated to be about 6.4 
1151 , i.e. thus about 17% larger than the associated experimental value of 5.5 [TTl . 

Nevertheless, and due to the potentially far-reaching consequences of the pro- 
vided new keV experimental results, we extended the aforementioned analysis 1131 
and Monte Carlo simulations 1 15 1 by including the strong dependence of the proton 
total cross section an on incident energy Eq f 14) into the calculations of double 
scattering. 

This dependence could be important, for the following reason. Neutrons which 
are not scattered once into the detector are multiply scattered. Such neutrons may 
have very high incident energies Eo, before they end up with final energy Ef = 
24.3 keV and become measured. Therefore they would have a much lower cross- 
section and hence the probability of scattering for these neutrons may be very 
low. Consequently it is conceivable that the multiple scattering contributions from 
protons could be significantly lower than those conventionally calculated on the 
basis of a fixed proton-neutron cross section (say, e.g., with cth = 15.7 bam, 
being valid for = 48.6 keV). 

Here it should be emphasized that the existing Monte Carlo routines for the 
calculation of multiple scattering effects cannot handle the present case of an 
depending on incident energy Eq; cf. 1 16|. The reason for this should be that, up 
to ca. 1 keV, the "free" total cross section of H remains constant, and the thus far 
existing NCS investigations were well within this energy range. 

To simplify the derivations and pinpoint their main physical aspects, we con- 
sider here double scattering from the water's protons only. Let Nh be the number 
of H atoms per unit volume and 



at "incident" neutron energy Ex before the a;-th neutron-proton collision. Let P be 
the scattering angle of the impinging neutron with initial energy Eq due to the first 
neutron-proton collision, and <j} being the associated azimuthal angle. (As usual, 
polar coordinates are used here.) Furthermore, let x be the angle between the two 
neutron-velocity vectors before and after the second scattering event. We consider 
neutrons with a total scattering angle 9 and with final energy Ef — 24.3 keV going 
into the detector. As is well known 1 17, 18|, the "intermediate" scattering angle x 
depends on both /3 and 6. A short calculation yields 



As conventionally, we consider here an infinite slab geometry for the sample, with 
thickness d, and the incident neutron beam in the direction z being perpendicular 
to the slab. 

According to standard theory f8l ll9l , the neutron energy Ei between first and 
second neutron-proton collision is 



fi{Ex) = NscTHiEx) 



(1) 



cos X = cos 6 cos /3 + sin 9 sin (3 cos (p 



(2) 



El = 



cos^ X 



(3) 



and the associated initial energy is 



Recall that here Ef = 24.3 keV. Thus, for some double-scattering events, Eo may 
be very large and the cross section gh{Eo) for the first scattering process very 
small fT4l . 

Furthermore we need the three attenuation factors (also called self shielding or 
absorption corrections) [20| for the neutron propagation through the sample: 

Uo = exp[-/i(£;o) za] (5) 

zo being the neutron path in the sample until its first collision, 

I7i =exp[-/i(£;i)si] (6) 

Si being the neutron path in the sample between its first and second collision, and 
finally 

-^^i^f) ZZTB 



Uf — exp 



(7) 



for the path between the second collision and the exit from the sample, in the fixed 
scattering direction 9 towards the detector. 

Furthermore, the probability for neutron scattering under scattering angle in 
the neutron-proton collision is proportional to cos V'; see Eq. (1.88) of 1191 . 

The intensity of impinging neutrons in the keV experiment is proportional to 

I{Eo) = So"° '' CD- 

As a result, the above partial processes lead straightforward to the following 
expression for the probability Pdse of a double-scattered neutron with final energy 
Ef = 24.3 keV measured with the detector at scattering angle 6: 

Pdse tx J I{Eo) [Uo fi{Eo) sin P cos /3] 

X [Ui n{Ei) cos x]Uf dzo dp d(j)dsi (8) 

(The squared brackets are for convenience only; they include the factors belong- 
ing to each scattering event.) The factor sin /? is due to the scattering after the 
first collision 1 19|. Note the presence of the factors /i(-Bo) and ^{Ei), which are 
associated with the scattering probability at the two scattering points. 



2.4 Main result - strongly anomalous R 

Note that all these equations are valid for the case of energy-dependent (jh{Eq). 
The probabilities for the twice scattered neutrons can be calculated with Eq. {Sj. 
The evaluation can be carried out with the Monte Carlo integration method; see 
e.g. Sect. 7.6 off2\\. 

The following observation is now of crucial importance. The same equation 
([8} can also be evaluated for the case affixed cross section oh =const. (For this, 
one simply has to fix the value of oh in Eq. 1.) Thus the results of both calcula- 
tions can be directly compared, as they are shown in Fig. 2. This comparison shows 
quantitatively the effect of the energy dependence of an, for various scattering an- 
gles. The results obtained with fixed cr// -value confirm the aforementioned Monte 
Carlo calculations of Mayers. The effect of the energy dependence aniEo) leads, 
at6 — 45°, to a slight decrease of the twice scattered neutrons of about 23 % only. 
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Figure 2: Intensity of double scattered neutrons from protons of pure H2O as a function 
of total scattering angle 9, calculated with Eq. (|8]l. Asterisks: Calculation with fixed 
cross section an =15.7 b. Circles: Energy-dependent (T^f(£'o), with Eq determined 
by Eq. (|4|i. Note the "forward" orientation of the twice scattered neutrons, which is 
equal in both calculations. The effect of the energy dependence of <th{Eo) leads to a 
decrease of the twice scattered neutrons of about 23 % at = 45°. This amounts to ca. 
3.5% of the complete intensity only; see the text. 



Furthermore, calculation of the single-scattered neutron intensity at was car- 
ried out by 

(9) 



(10) 



Psseoc / I{Eo)[U;,fiiE'o)cose]U'fdzo 
Jo 

with E'o = Ef/cos^ 9, Uq = exp[-^(£'o) zq] and 



U'f = exp 



(This integral can be calculated analytically, too 1201 .') Comparison of the above 
results showed that the ratio of double- to single-scattered neutrons is about 0.15 
at S = 45°. Thus we may proceed to the following crucial conclusion: If one 
approximately neglects higher-order scattering events (which were estimated to 
be about 2-3% of the total scattering), then the aforementioned 23%-reduction 
of the double-scattered neutrons amounts to ca. 3.5% of the complete intensity at 
B = 45° . This furthermore implies that - neglecting double-scattering from D2O - 
the earlier calculated ratio J? = 6.78 ( 13| should be slightly reduced to R = 6.55, 
which is still about 19%, and thus significantly, larger than the experimentally 
measured ratio Rexp = 5.5 

Moreover, if one would also take into account double scattering events from 
D2O, it is obvious that the considered ratio R of integrated scattering intensities 
from H2O and D2O will be slightly increased again. Thus we may safely conclude 
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that the theoretically predicted value of the ratio R still remains about 20% larger 
than the experimental one Rexp = 5.5. 

This crucial result is fully in line with the previous calculations based on single 
scattering events only 1 13 1. It shows that inclusion of double scattering, as well as 
the strong dependence of the proton total cross section an on incident energy Eo 
into the calculations, confirm again the anomaly of the intensity ratio R. As the 
overall errors of the keV-measurements were about 3% while the actual statistical 
errors were about 2% [llj . one concludes that the revealed anomaly of ca. 20% is 
highly significant. 

3 Theoretical remarks 

Originally it has been proposed 1 1 1 that the considered effect is caused by short- 
lived and spatially restricted entanglement. Since the typical neutron- and electron- 
proton interaction time in NCS and ECS (i.e. Tsc) lies within the attosecond range, 
it is expected that decoherence may still not be fully effective. Published theoreti- 
cal models attribute this effect to: 

[A) Modification of scattering due to "identity of particles" in the scattering 
system. The influence on scattering of entanglement of the spin and spatial degrees 
of freedom of identical particles (i.e., quantum exchange correlations) has been 
stressed in Refs. [22 1. In particular, NCS from pairs of protons and deuterons has 
been calculated. 

(B) Contribution of electronic degrees of freedom to the dynamics of a struck 
proton (deuteron) interacting with its environment. E.g., breakdown of the Bom- 
Oppenheimer (BO) approximation in the final state of the NCS process and (-B.1) 
additional excitations of the electronic system (23l,'24|, and/or (-B.2) decoherence 
accompanying short-lived spatial entanglement of a struck proton with adjacent 
electrons and perhaps also nuclei |25|. In the models of this category, spin entan- 
glement and/or quantum exchange correlations between identical particles play no 
role. 

To test certain contradictory predictions of these theories, we recently mea- 
sured by NCS (a) the equimolar H2-D2 mixture and (6) the mixed-isotope sys- 
tem HD (liquids, both at 20 K). The crucial result was that both systems reveal 
the same, strong anomalous shortfall (about 30%) of the ratio an / of H and 
D cross-sections 1261 . Since HD exhibits no exchange correlations, this result 
demonstrates for the first time that these correlations play no significant role in this 
effect, thus refuting corresponding theoretical models claiming its interpretation. 
This conclusion is in line with a recently presented full calculation of the scatter- 
ing function S{q, ui) by Sugimoto et al. [27J , who found that indistinguishability 
of particles cannot represent the physical origin of the observed effect. 

In contrast, our findings 1 26 1 are consistent with alternative theoretical models 
of category (B), in which attosecond dynamics of electronic degrees of freedom 
— e.g. via violation of the BO approximation — is considered to participate sig- 
nificantly to the dynamics of an elementary neutron-proton (-deuteron) scattering 
process. The scattering times of NCS and ECS are similar to the characteristic 
time of "electron motion", so that the Bom-Oppenheimer approximation is not ap- 
plicable here; cf. Ref. [2a]. However, it should be stressed that a detailed treatment 
of the BO failure and a quantitative estimate of its possible contribution to the con- 
sidered NCS-anomaly does not exist yet, as convincingly discussed in Ref. (28). 
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4 Additional remarks 



The preceding derivations concerning the new keV experimental results II II have 
various far-reaching theoretical and experimental consequences. Some of them are 
as follows. 

First, it should be stressed that the intensity ratios determined in this experi- 
ment can only provide information about a possible difference between the scat- 
tering behavior of H and D, and not about that of H alone. In view of several 
theoretical works P 23| 1241 |25l , it may be expected that, due to the high energy 
transfers, both protons and deuterons should violate the BO approximation, thus 
leading to the above "anomaly" for both H and D. Thus, the measurement of the 
scattering intensity ratio Rexp of H2O relative to that of D2O is not appropriate 
for the search of the anomalous scattering effect under consideration. Obviously, 
an amended keV experiment would be to measure the scattering intensity ratio of 
H2O (D2O, and their mixtures) relative to that of a heavy metal sample, e.g. Pb. 
This experiment was recently proposed by the author. 

Second, already earlier NCS from D-containing materials have shown that 
also D exhibits a small anomalous shortfall of scattering intensity. E.g. in [2c] 
was reported that this shortfall in NbDo.8 was about 10%. Furthermore, several 
NCS measurements on pure D2O showed up a shortfall of the cross sections ra- 
tio aojao of about 10-15% (unpublished data). Thus, the aforementioned 20%- 
anomaly of the measured ratio i?ea;p in the keV experiment indicates that the inte- 
grated scattering intensity of pure H2O may exhibit an anomalous shortfall larger 
than 20%. In view of these considerations, the aforementioned proposed experi- 
ment may be expected to yield an anomalous intensity ratio of H2O to Pb of the 
order of 25% and more. 

Summarizing, we conclude that the considered scattering effect is present at 
both 5-100 eV 1 1| and 24-150 keV 1 11 1 ranges of incident energies. Obviously, 
the novel experiment 1 1 1 1 and its correct analysis established above open up new 
perspectives for neutron research on the above attosecond effect, and thus they may 
have far reaching consequences for current and future experimental and theoretical 
investigations. 
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